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A condition on the Hamiltonian of a time-dependent quantum mechanical system is derived which, 
if satisfied, impHes optimal adiabaticity (defined below). The condition is expressed in terms of the 
Hamiltonian and in terms of the evolution operator related to it. Since the latter depends in a 
complicated way on the Hamiltonian, it is not yet clear how the condition can be used to extract 
useful information about the optimal Hamiltonian. The condition is tested on an exactly-soluble 
time-dependent problem (a spin in a magnetic field), where perfectly adiabatic evolution can be 
easily identified. 



I. INTRODUCTION 

The adiabatic theorem in quantum mechanics was de- 
veloped when quantum mechanics was still in its infancy 
[1]. (See p] for an early English-language reference, and 
[3] for a recent overview.) It is used in virtually every 
area of quantum physics. The essential idea underly- 
ing the theorem is very simple. Let the Hamiltonian 
be H{t) and define its instantaneous eigenstates |n(i)) 
and eigenenergies En{t) as the solutions to the time- 
independent Schroedinger equation H(t)\'ip) ~ E\ijj) 
(where t is viewed as a parameter). The states \n{t)) 
evolve in time, as does the solution \^p{t)) of the time- 
dependent Schroedinger equation 

\m) = m \m) ■ 

In general there is no obvious relation between these 
two evolutions. The adiabatic theorem states that the 
system follows the instantaneous eigenstates for infinitely 
slow evolution. More precisely, suppose the system starts 
in |n(0)) at t = 0, and that \n{t)) is gapped at all times 
(that is, \En{t) - E„,{t)\ > Eq > for all t and for 
all m ^ n). Suppose that we can control the speed with 
which the Hamiltonian varies by making the replacement 
H{t) — H\{t) = H{Xt) with the time interval changed 
accordingly, [0,T] [0,T/A], so that the limit A ^ 
corresponds to infinitely slow evolution. Then the adia- 
batic theorem states that the solution to the Schroedinger 
equation 

if^\Mt)) = H,it)\Mt)), \MO)) = \nm 

satisfies 

lim \^x{T/\)) = (phase) |n(T)) + 0(A). 

A-)-0 

Thus, as A — > 0, the probability for the system to make 
a transition to a different instantaneous eigenstate goes 
to zero. We then say that the evolution is adiabatic. 



Infinitely slow evolution is time consuming, to say the 
least, so it is of interest to be able to make a statement 
about finitely slow evolution. Intuitively, if the evolution 
is sufficiently slow, the adiabatic approximation 

- (phase) \n{t)} (1) 

should be reasonable. 

This begs the question: how slow must the evolution 
be for ([T]) to be a good approximation? In a situation 
where adiabatic evolution is desirable, it is obviously im- 
portant to know just how slowly the system must evolve 
in order for this to be the case, in order to accomplish the 
task required as quickly as possible. (In the case of adi- 
abatic quantum computing [HE], for instance, how the 
slowness required scales with the size of the system being 
studied enables comparison of this approach to conven- 
tional circuit-based quantum computing.) 

There are two time scales at play: the time scale of 
the evolution of the Hamiltonian and a scale related to 
the difference in energies. The former is often estimated 
to be of order | {m{t)\ h{t))\~^; the latter can be written 
\Ejn — En\~^. Usually, the evolution can be considered 
slow [5] if 



E — E 

But this condition does not guarantee that ([l]) is a good 
approximation. First, it does not necessarily imply slow 
evolution [7j. Secondly, even in circumstances where it 
does imply slow evolution, one cannot conclude that the 
adiabatic approximation is true indefinitely; rather, it 
implies that the system will "escape" from \n{t)) more 
slowly than for a "generic" (non-slow) evolution. Much 
work has been done recently studying the adiabatic ap- 
proximation and attempting to give one or more condi- 
tions guaranteeing its validity [71-1^ . 

In this paper, our goal will be to devise a criterion 
which implies that a Hamiltonian is optimal in the sense 
that the evolution is optimally adiabatic. If \ip{0)) = 



2 



|n(0)), then we define the adiabaticity as 
A{t)^\{nit)\m)f- 



(2) 



A cannot exceed unity, and A — 1 corresponds to per- 
fectly adiabatic evolution (that is, evolution where the 
adiabatic approximation is exactly satisfied). 

The specific question we will address is as follows. Sup- 
pose we are given a physical system described by a time- 
dependent Hamiltonian, and that the duration of evo- 
lution (call it T) and the initial and final Hamiltonians 
{H{0) and H{T)) are given. Suppose furthermore that, 
as above, the system is initially in the state |n(0)). The 
question is: Can we determine an interpolating Hamilto- 
nian H(t) which maximizes the (final) adiabaticity, 



A[H{t)] = \{n{T)\i;{T))Y 



(3) 



In what follows, we will devise such a criterion, and will 
test it on an exactly-soluble system. Unfortunately, as 
will be seen below, the criterion is expressed in terms of 
the evolution operator. Thus, it is not clear how it can 
be used on a system for which this operator cannot be 
determined. In the next section we begin by observing 
that there is a trivial, but impractical, way of obtaining 
perfectly adiabatic evolution. We then restrict and for- 
mulate the question addressed above. In the following 
section we explain the optimization procedure which re- 
sults in out main result, ([6|. Following this, we test the 
result on what is probably the simplest case of a time- 
dependent Hamiltonian, one which can be solved exactly: 
a spin in a rotation magnetic field. We end with a sum- 
mary and concluding remarks. 



II. STATEMENT OF THE PROBLEM 



so in this limit the state of the system at the end of this 
stage is a phase times \n{T)). During the third step, the 
Hamiltonian is again multiplied by a linear function of 
time, going from KHi to Hi, and as in the first stage the 
state of the system only changes by a phase. 

Suppose we restrict ourselves to situations where the 
energy eigenvalues are constant, -E„(t) = i?„; this elim- 
inates the trivial way just described, and simplifies the 
analysis to follow. The Hamiltonian then evolves via a 
unitary transformation 

H{t)^V{t)H{Q)V\t), 

where V^{t)V{t) — 1. We can assume V{Q) = 1, and the 
final value V{T) is fixed since the final Hamiltonian is 
presumed to have been specified in advance; the instan- 
taneous eigenstates are \n{t)) = V{t)\n{0)). We write 
A[Hit)] ^ A[Vit)]. 

Let the time evolution operator associated with H{t) 
be U{t\ so the state is given by |i/i(t)) = U{t)\n{0)). 
Then (Isl) becomes 



A[Vit)] = \{nmvHT)U{T)\nm\ 



(4) 



and we wish to find a condition on V{t) for which the 
adiabaticity is maximal. 



III. OPTIMIZATION 

We will adopt a variational approach to find a condi- 
tion on V{t). Suppose that the adiabaticity is optimized 
for a certain matrix Vo{t) (to be determined). Then for 
any variation of V{t) about Vo(i) which is zero initially 
and finally (so that the initial and final Hamiltonians are 
unaffected), A is stationary to first order: 



Before we begin, note that there is actually a trivial 
(although impractical) way of attaining perfectly adia- 
batic evolution Suppose Hg (where s : — >■ 1) is 
a sufficiently smooth family of Hamiltonians which in- 
terpolates between H{0) and H{T) (so that Hq = H{0) 
and Hi — H{T)) and for which the instantaneous eigen- 
state jris) is always gapped. Then the following time- 
dependent Hamiltonian will give perfectly adiabatic evo- 
lution in the limit A — >■ cx): 



HK[t) 



3A- 2 



l + f(A-l))i/o t 



'-l + f)AHi t 
f(l-A))i/i t 



- 
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h 
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The evolution is divided into three steps. During the first 
step, the Hamiltonian is simply multiplied by a linear 
function of time, going from Hq to AiJp; |?^(0)) remains 
an eigenstate of the Hamiltonian so the state of the sys- 
tem only changes by a phase. During the second step, the 
Hamiltonian evolves (again linearly) from KHq to Ai?i; 
as A goes towards infinity, the probability of making a 
transition from \n{t)) to a different state drops to zero. 



5A 



5V{t) 



0. 



V{t) = Vo{t) 



(Of course, solutions to this equation will be local ex- 
trema, not necessarily global minima, of the adiabaticity.) 

We write Hoit) = Vo{t)H{0)Vi^ (t), with Uo{t) the cor- 
responding evolution operator and Aq the optimal (as- 
sumed maximal) adiabaticity. 

Now consider a small variation of V{t): 

Vit) = {l + th{t))Vo{t), 

where h{t) = h^{t) (in order for V{t) to be unitary), 
h{t) <^ 1 (meaning the matrix elements are much less 
than 1) and h{0) = h{T) = 0. Then 



SA 



SVit) 
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V{t)=V„{t) 



5h{t) 



h{t)=0 



In Q, V{T) = Vo{T) is independent of h{t) (so in what 
follows we will write (n(0)|y^(T) ^ {n{T)\), but U{T) 
depends non-trivially on h{t). 
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The change Vo(i) — ^ ^(i) will induce a small change 
in the evolution, so we can write 

U{t) = (l + ^A:(^))C/o(^), 

where k{t) = kHt), fc(0) = and A: < 1. Then 

A^\{n{T)\{l + ik{T))Uo{T)\n{0)) 

= + 2ilm {{n{T)\ k{T)Ua{T) |n(0)) x 

{nmuUT)\niT))]+0{k^). 

That yl is stationary implies 

Im{(n(r)| k{T)UoiT) |n(0)) (n(0)| U^T) |n(r))} = 

(5) 

for any h{t). 

We must now express k{T) in terms of h{t). We note 
that for the unperturbed problem 

ij^Uoit) = Ho{t)Uo{t), Uo{0) = 1, 
while for the perturbed problem 

ijU{t)=H{t)U{t), C/(0) = 1, 

with the Hamiltonians and evolution operators related 
by 

H{t) = (1 + ih{t))Ho{t){l - ih(t)), 
U{t) = {l + ikit))Uo{t). 

Direct substitution and taking linear terms in both h and 
k yields 



[Ho,k]~[Ho,hl fc(0) = 0. 



Standard techniques yield the following solution, as can 
be verified directly: 



kiT)^Uo{T)i / dtU'oit)[Hoit),hit)]Uom'o{T) 
Jo 

Substituting in ([s]), we see that 

Re dt { (n(T) | C/q (T) Ul (<) [H^ {t) , h{t)] C/o (i) |n(0)) 



X {nmul{T)\n{T))] 



= 



for any Hermitian h{t). Each of the factors in the braces 
can be analyzed semi-intuitively. The second, which can 
be loosely described as the "square root" of Aq , is simply 
the amplitude for the optimal evolution to be adiabatic. 
The first factor, read right to left, is: starting in the 
initial state, evolving optimally for a time an effect of 
the perturbation acting on the state at time i, evolving 
for the remaining time T — t, and projecting onto the 
final instantaneous eigenstate. Thus, roughly speaking. 



the sum of all possible first-order changes to A must be 
zero. 

We can progress further by analyzing what we mean 
by the statement "for any Hermitian h{t)." We can write 
h{t) = Xifi{t) (sum on i implicit), where {A^} are a ba- 
sis of Hermitian matrices of the appropriate dimension 
and fi{t) are arbitrary functions (except that they are 
zero initially and finally). We can vary these functions 
independently, and indeed if one of them is nonzero at 
one instant and the rest are zero for all times, the above 
condition becomes 



Re 



{(n(T)l UoiT)U^oit)[Hoit),X,]Uo{t) |n(0)) 

x(n(0)|C/o^(T)Kr))}=0, (6) 

an equation which must be true for all intermediate times 
t and for all i. This is in principle an equation to be 
solved for Hq. However, as mentioned above, Uq depends 
in a highly non-trivial way on Hq. Thus it is far from 
obvious how to use it to learn something about a system 
for which the time evolution operator is unknown (and 
if it is known, then the adiabatic approximation is of 
limited use). 

IV. TEST OF THE CONDITION 

We can at least test the condition on an exactly soluble 
model: the much-studied case of a spin-1/2 particle in 
a constantly rotating magnetic field. The geometry is 
illustrated in Figure [TJ 




FIG. 1. Rotating magnetic field with angular frequency uj and 
angle with respect to axis of rotation 9. 

The Hamiltonian is 

UJQ 



(7) 



where n is a unit vector in the direction of the magnetic 
field and is the Larmor frequency. The time evolution 
operator is 



(8) 
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where /3, w are defined in Figure [2] 

Since the model can be solved exactly, the adiabaticity 
can be calculated analytically as a function oft; the result 
(supposing we have started in an eigenstate of the initial 
Hamiltonian, e.g. exp(— i6'cr2/2) |+)) is 

A is obviously maximized when ujt = 27rm, where m G Z. 

Thus, we can imagine being asked to find a Hamilto- 
nian which optimizes the adiabaticity for the following 
parameters of the problem: 

• system: spin 1/2 

• H{T) = -5^iie~*'#'"3e"*5'^^cr3e*2'"2e''^'^3 

• T — 2nnr/u) with w as defined in Figure [2] 

In principle, we would like to use ^ to determine an 
optimal Hamiltonian Hq (t) . We have not found a way to 
do so directly, but we can at least check that ([t]) is indeed 
optimal. This is simply a matter of substituting ([t]) and 
([s]) into (|6|. This is somewhat tedious but absolutely 
straightforward, and we see that indeed for T = 2m7r/aj, 
([t]) does indeed result in a solution of Q, as expected. 



V. CONCLUSIONS 

A condition was derived which implies that the adia- 
baticity (defined in ([2|) is maximal. This condition in 
principle determines which Hamiltonian or Hamiltonians 
give rise to perfectly adiabatic evolution. However the 
condition is expressed in terms of the evolution operator, 
which depends on a complicated way on the Hamilto- 
nian, so it is not yet clear how to extract useful infor- 
mation about a system for which the evolution operator 
cannot be determined. We verified that the optimization 
condition is indeed satisfied when the evolution is opti- 
mal, in the case of an exactly-solvable system: a spin in 
a uniformly rotating magnetic field. 

We thank Manu Paranjape for useful discussions. This 
work was funded in part by the Natural Science and En- 
gineering Research Council of Canada. 
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FIG. 2. Illustration of /?, 
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